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LOCAL TO GLOBAL ALGORITHMS FOR THE GORENSTEIN 
ADJOINT IDEAL OF A CURVE 


JANKO BOHM, WOLFRAM DECKER, SANTIAGO LAPLAGNE, AND GERHARD PFISTER 


Abstract. We present new algorithms for computing adjoint ideals of curves 
and thus, in the planar case, adjoint curves. With regard to terminology, we 
follow Gorenstein who states the adjoint condition in terms of conductors. 

Our main algorithm yields the Gorenstein adjoint ideal 0 of a given curve as 
the intersection of what we call local Gorenstein adjoint ideals. Since the respec¬ 
tive local computations do not depend on each other, our approach is inherently 
parallel. 

Over the rationals, further parallelization is achieved by a modular version of 
the algorithm which first computes a number of the characteristic p counterparts 
of 0 and then lifts these to characteristic zero. As a key ingredient, we establish 
an efficient criterion to verify the correctness of the lift. 

Well-known applications are the computation of Riemann-Roch spaces, the 
construction of points in moduli spaces, and the parametrization of rational 
curves. 

We have implemented different variants of our algorithms together with Mnuk’s 


approach Mnuk 1997 in the computer algebra system Singular and give timings 


to compare the performance of the algorithms. 


1. Introduction 


In classical algebraic geometry, starting from Riemann’s paper on abelian func¬ 
tions [Riem ann 1857], the adjoint curves of an irreducible plane curve T have been 
used as an essential tool in the study of the geometry of T. The defining property 
of an adjoint curve is that it passes with “sufficiently high” multiplicity through the 
singularities of T. There are several ways of making this precise, developed in clas- 


and in more recent work by [Grbbner 1941 

Gorenstein 1952 and [van der Waerden 

1939, 

Keller 1974 . We refer to [Keller 1965 , Greco and Valabrega 1979 , |Greco 

and Valabrega 1982 , and Ciliberto and Orecchia 1984 for results comparing the 


different notions: whereas the adjoint condition given by Brill and Noether is more 
restrictive, the notions of adjoint curves given by the other authors above coincide. 

In this paper, we always consider adjoint curves in the less restrictive sense. In 
fact, we rely on Gorenstein’s algebraic definition which states the adjoint condition 
at a singular point P E T by considering the conductor of the local ring 0r,p in its 
normalization. It is a well-known consequence of Max Noether’s Fundamentalsatz 
that the adjoint curves of any given degree m cut out, residual to a fixed divisor 
supported on the singular locus of T, a complete linear series. Of fundamental 
importance is the case m = degT — 3 which, as shown by Gorenstein, yields the 
canonical series. 

The ideal generated by the defining forms of the adjoint curves of T is called 
the adjoint ideal of T. In |Arbarello and Ciliberto 1983|, the concept of adjoint 
ideals is extended to the non-planar case: consider a non-degenerate irreducible 
curve r CPJ = Proj(S'), and let I be a saturated homogeneous ideal of S which 
is supported on the singular locus of T. Then, roughly speaking, / is an adjoint 
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ideal of F if its homogeneous elements of degree m 3> 0 cut out, residual to a 
fixed divisor supported on the singular locus, a complete linear series. As pointed 
out in |Arbarello and Ciliberto 1983j, the existence of adjoint ideals is implicit in 
classical papers: examples are the Castelnuovo adjoint ideal and the Petri adjoint 
ideal. In [ Ciliberto and Orecchia 1984], it is shown that Gorenstein’s condition leads 
to the largest possible adjoint ideal, containing all other adjoint ideals, and now 


referred to as the Gorenstein adjoint ideal © = ©(r). See Ciliberto and Orecchia 


1984 for some remarks on how the different concepts of adjoint ideals compare in 


the non-planar case. 

With regard to practical applications, adjoint curves enter center stage in the clas¬ 
sical Brill-Noether algorithm for computing Riemann-Roch spaces, which in turn can 
be used to construct Goppa codes (see |Le Brigand and Risler 1988) ). Furthermore, 
linear series cut out by adjoint curves allow us to construct explicit examples of 
smooth curves via singular plane models; a typical application is the experimental 
study of moduli spaces of curves. If the geometric genus of a plane curve T is zero, 
then the adjoint curves of degree degT — 2 specify a birational map to a rational 
normal curve. Based on this, we can find an explicit parametrization of T over 
its field of definition, starting either from the projective line or a conic. See |Bohm 


1999, Bohm et al. 2015c and the implementation in the Singular library Bohm et 


al. 2012c . Algorithms for parametrization, in turn, have applications in computer 


aided design, for example, to compute intersections of curves with other algebraic 
varieties. See also [Sendra et al. 


2008 


A well-known algorithm for computing the Gorenstein adjoint ideal © = ©(r) in 

This algorithm makes use of linear algebra to 


the planar case is due to Mnuk 1997 


obtain © from an integral basis for the normalization k[C\, where C is an affine part 
of r containing all singularities of T. Efficient ways of finding integral bases rely on 
Puiseux series techniques (see [van Hoeij 1994 , Bohm et al. 2015a|). This somewhat 


limits Mnuk’s approach to characteristic zero. The same applies to the algorithm 

which also computes the Gorenstein adjo int i deal 


of El Kahoui and Moussa 2014 


of a plane curve from an integral basis of k[C\. The approach of Orecchia and 
Ramella, 2014|, on the other hand, is limited to ordinary multiple points. 


In this paper, we present a new algorithm for computing ©. This algorithm is 
highly efficient and not restricted to the planar case, special types of singularities or 
to characteristic zero. The basic idea is to compute © as the intersection of “local 
Gorenstein ideals”, one for each singular point of T. Each local ideal is obtained 
via Grobner bases, starting from a “local contribution” to the normalization k[C] 
at the respective singular point. To find these contributions, we use the algorithm 
from | Bohm et al. 2012a| which is a local variant of the normalization algorithm 


designed in Greuel et al. 2010a . In practical terms, given any field of definition 


L C k, we treat the points in a complete set of conjugate singularities simultaneously. 

Our approach is already faster per se. In addition, it can take advantage of 
handling special classes of singularities in an ad hoc way. Above all, it is inher¬ 
ently parallel. For input over the rationals, further parallelization is achieved by a 
modular version of our algorithm which first computes a number of characteristic 
p counterparts of © and then lifts these to characteristic zero. To apply the gen¬ 
eral rational reconstruction scheme from |Bohm et al. 2012dj, we prove an efficient 
criterion to verify the correctness of the lift. 

Our paper is organized as follows: In Section [2j we discuss algorithmic normaliza¬ 
tion. In Section [3j we review the definition of adjoint ideals and some related facts. 
In Section [4j we describe global algorithmic approaches to obtain ©. We first discuss 
Mnuk’s approach. Then we describe a global approach which relies on normalization 
and Grobner bases. In Sections [5] and [6| we present our local to global algorithm for 
finding © via normalization and Grobner bases. Section [7] pays particular attention 
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to the planar case, commenting on the direct treatment of special types of singu¬ 
larities. In Section [8j we discuss the modular version of our algorithm. Finally, in 
Section [9j we compare the performance of the different approaches, relying on our 
implementations in the computer algebra system SINGULAR, and running various 
examples coming from algebraic geometry. 


2. Algorithms for Normalization 


We begin with some general remarks on normalization and the role played by the 
conductor. For these, let A be any reduced Noetherian ring, and let Q(A) be its 
total ring of fractions. Then Q(A) is again a reduced Noetherian ring. We write 

Spec(A) = {P C A | P prime ideal} 

for the spectrum of A. The vanishing locus of an ideal J of A is the set V(J) = 
{P G Spec(A) | P d J}. 

The normalization of A, written A, is the integral closure of A in Q(A). We call 
A normalization-finite if A is a finite A-module, and we call A normal if A = A. 
We denote by 

N(A) = {P G Spec(A) | Ap is not normal} 
the non-normal locus of A, and by 

Sing(A) = {P G Spec(A) | Ap is not regular} 
the singular locus of A. 


Remark 2.1. Note that N(A) C Sing(A). Equality holds if A is of pure dimension 
one. Indeed, a Noetherian local ring of dimension one is normal iff it is regular 
(see (de Jong and Pfister 2000, Thm. 4.4.9]). 

Definition 2.2. If R C S is an extension of rings, the conductor of A in B is 

Cs/R = {r e R\rS C S} . 

Note that Cp/p is the largest ideal of R which is also an ideal of S. 

Notation 2.3. If A is a reduced Noetherian ring as above, we write 

C A = Ca/a = {« £ A | aA c A}. 

Lemma 2.4. We have N(A) C V{Ca)- Furthermore, A is normalization-finite iff 
Ca contains a nonzerodivisor of A. In this case, N(A) = V(Ca)- 

Proof. See [Greuel and Pfister 2008, Lemmas 3.6.1, 3.6.3]. □ 

Remark 2.5 (Splitting of Normalization). Finding the normalization can be reduced 
to the case of integral domains: If P\ ..., P s are the minimal primes of A, then 


A ^ A/Pi x ■ • • x A/P s 
(see jde Jong and Pfister 2000, Thm. 1.5.20]). 

Remark 2.6. Let k be a field. An affine k-domain is a finitely generated fc-algebra 


which is an integral domain. By Emmy Noether’s finiteness theorem (see Eisenbud 


1995, Cor. 13.13]), any such domain is normalization-finite, and its normalization is 


an affine /c-domain as well. Geometrically, by gluing, this implies that any integral 
variety X over k admits a (unique) normalization map X —> X, where X is again 
an integ ral v ariety over k (see, for example, Liu 2002, Sec. 4.1.2]). Specifically, by 
Remark 


2.1 


if r is a curve over k, we get the nonsingular model n : T —> T. 
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Now, we briefly discuss algorithmic normalization, 
normalization algorithm of Greuel, Laplagne, and Seelisch 
which is an improvement of de Jong’s algorithm (see [de Jong 1998 ~ [Decker et al. 
19991). This algorithm, to which we refer as the GLS Algorithm, is based on the 


We begin by recalling the 
Greuel et al. 2010a], 


normality criterion of Grauert and Remmert. To state this criterion, we need: 


Lemma 2.7. Let A be a reduced Noetherian ring, and let J C A be an ideal which 
contains a nonzerodivisor g of A. Then: 

(1) If p 6 Hom^J, J), the fraction p{g)/g £ A is independent of the choice of 
g, and p is multiplication by p(g)/g. 

(2) There are natural inclusions of rings 

A C HomAJ, J) = ~(gJ ’-a J) T A C Q(A), a i-A p a , P ^ , 

9 g 

where p a is multiplication by a. 


Proof. See |Greuel and Pfister 2008 


Lemmas 3.6.1, 3.6.3]. 


□ 


Proposition 2.8 (Grauert and Remmert Criterion). Let A be a reduced Noetherian 
ring, and let J C A be a radical ideal which contains a nonzerodivisor g of A and 
satisfies V(Ca) C V(J). Then A is normal iff A = Hom^J, J) via the map which 
sends a to multiplication by a. 


Proof. See Grauert and Remmert 1971 , [Greuel and Pfister 2008, Prop. 3.6.5]. □ 


Definition 2.9. A pair (J, g) as in the proposition is called a test pair for A, and 
J is called a test ideal for A. 


If k is a field and A is an affine &:-domain, then test pairs exist by Lemma 2.4 and 
Emmy Noether’s finiteness theorem. If, in addition, k is perfect, a test pair can be 
found by applying the Jacobian criterion (see [Eise nbud 1995, Thm. 16.19] for this 
criterion). In fact, in this case, we may choose the radical of the Jacobian ideal M 
together with any nonzero element g of M as a test pair. Given a test pair ( J,g ), 
the basic idea of finding A is to enlarge A by a sequence of finite extensions of affine 
A;-domains 

A i+ 1 = Horn Ai {Ju Ji ) = ~{gJi -A, Ji)cAc Q(A), 

with Aq = A and Ji = \J J A 
stop. According to 

bUi 


until the Grauert and Remmert criterion allows one to 


Greuel et al. 2010 a 


each Ai can be represented as a quotient 
C Q(A), where Ui C A is an ideal and di £ Ui is nonzero. In this way, all 


computations except those of the radicals Ji may be carried through in A. 


Example 2.10. For 

A = C[x, y} = C[X, Y}/(X 5 - Y 2 (Y - 1 ) 3 >, 
the radical of the Jacobian ideal is 


J ■= (x,y(y - i))^ > 

so that we can take (J, x) as a test pair. Then, in its first step, the normalization 
algorithm yields 

Ai = V, =±(x,y(y- l) 2 ) A . 

In the next steps, we get 

a 2 = ^ u 2 = ~2 ( x 2 > x v(y ~ VMy - !) 2 )a 

and 

^3 = ^ 3 U 3 = ^ (x 3 ,x 2 y{y - 1 ),xy(y - l) 2 , y 2 (y - l) 2 ) A . 

In the final step, we find that A 3 is normal and, hence, equal to A. 
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Next, we describe a local to global variant of the GLS algorithm, given in |Bohm et 


al. 2012a|, which is a considerable enhancement of the algorithm, and which serves 


as a motivation for our local to global approach to compute the Gorenstein adjoint 
ideal . This variant is based on the following two observations from |Bohm et al. 


2012a : First, the normalization A can be computed as the sum of local contributions 
a« C A, and second, local contributions can be obtained efficiently by a local 


Ac 


variant of the GLS algorithm. For our purposes here, it is enough to present the 
relevant results in a special case. Here, as usual, if P is a prime of a ring R, and M 
is an R- module, we write Mp for the localization of M at R \ P. 

Proposition 2.11. Let A be an affine k-domain of dimension one, and let Sing(A) = 
{Pi,..., P s } be its singular locus. For i = 1,... ,s, let an intermediate ring A C 
yjb) c A be given such that A p = Ap i . Then 

J2 A (i) = A. 

i= 1 

Proof. See [Bohm et al. 2012a, Prop. 15]. □ 

Definition 2.12. A ring A W as above is called a local contribution to A at P % . It 
is called a minimal local contribution if A p = Ap j for j R i. 

The computation of local contributions is based on the modified version of the 
Grauert and Remmert criterion below: 

Proposition 2.13. Let A be an affine k-domain of dimension one, let A C A' be a 
finite ring extension, let P G Sing (A), and let J' = \/ PA!. If 

A' 9* Horn A '{J',J') 

via the map which sends a' to multiplication by a!, then A' p is normal. 

Proof. See (Bohm et al. 2012a, Prop. 16]. □ 

Considering an affine domain A of dimension one over a perfect field k , let P G 
Sing(A). Choose P together with a nonzero element g in P instead of a test pair as 
in Definition 2.9 Then, proceeding as before, we get a chain of affine fc-domains 

A C Ai C • • • C A m C A 
such that A m is a local contribution to A at P. 

Remark 2.14. Given A as above, a finite ring extension A C A', and a prime P € 
Sing(A), let Q € Sing(A) be a prime different from P, and let J' = \JPA'. Then 

HomA'(J', J')q = Horn a’ q (Jq, Jq) 

— Hom a' q (Aq, Aq) = A'q 

(see |Eise nbud 1995, Proposition 2.10]). Inductively, this shows that the algorithm 
outlined above computes a minimal local contribution to A at P. Note that such 
a contribution is uniquely determined since, by definition, its localization at each 
P G Spec(A) is determined. 

Example 2.15. In the case of Example |2.10 there are two singularities P\ = 
(x,y) and P 2 = {x,y — 1). For Pi, the local normalization algorithm yields Ap 1 = 
(^CA)pi, where 

di = x 2 and IR = (x 2 , y(y-lf) A . 

For P- 2 , we get Ap 2 = (±U 2 )p 2 , where 

d 2 = x 3 and U 2 = (x 3 , x 2 y 2 (y - 1), y 2 {y - l) 2 ^ . 
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Combining the local contributions, we get 


-U = + ^U 2 

a cl\ d 2 


with d = x 5 and 


U = (x 3 , xy(y - l) 3 , x 2 y 2 (y - 1), y 2 (y - l) 2 ^ . 


Note that U coincides with the ideal U 3 computed in Example 2.10 


Notation 2.16. In our applications, A will always be the coordinate ring k[C} = 
k [X\, ...,X r \ /1(C) of an integral affine curve C C A r k over a perfect field k. Given 
a poin10 P € C, by abuse of notation, if I C k[xi,...,x r \ is an ideal properly 
containing 1(C), we will write Ip for the ideal of the local ring Oc,p obtained by 
mapping I to k[C] and localizing at P. Likewise for the homogeneous localization 
of a homogeneous ideal in the projective case. 


3. Adjoint ideals 


Let k be a field, and let T C V r k be an integral non-degenerate projective curve. 
Write S = k [Aq,..., X r ] for the homogeneous coordinate ring of P£, I(r) C S for the 
homogeneous ideal of T, fc[r] = S/I(T) for the homogeneous coordinate ring of T, 
and Sing(r) for the singular locus of T. 

Let 7r : r — y r be the normalization map, let P be a point of T, and let 0p,p be 
the local ring of T at P. Then the normalization Op,p is a semi-local ring whose 
maximal ideals correspond to the points of T lying over P. Furthermore, Op,p is 
finite over Op,p and, thus, a finite-dimensional fc-vector space. The dimension 

<5p(r) = 5((Dr,p) = dim*. Or,p/Or,p 

is called the delta invariant of T at P. The arithmetic genus of T is p a (T) = 1 —Pr(0), 
where Pr is the Hilbert polynomial of k[T]. Making use of the (global) delta invariant 

^(r) = £p e sin g( r)Mr) 
of T, the geometric genus p(T) of T is given by 


(see 


Hironaka 1957 


p(T) = p(T) = p a (T) ~ S(T) 

). If T is a plane curve of degree n, we have p a (r) = (” 2 1 )- 


Following the presentation in Chiarli 1984], we now recall the definition and char¬ 


acterization of adjoint ideals due to [Arbarello and Ciliberto 1983] and [Ciliberto and 


Orecchia 1984 


Let / = © m>0 dm C S = k[X 0 ,..., X r \ be a saturated homogeneous 
ideal properly containing /(T). Pulling back Proj(S'/J) via n, we get an effective 
divisor A (I) on T. Let H be a divisor on T given as the pullback of a hyperplane 
in Pf. Then, since any divisor on T cut out by a homogeneous polynomial in I is of 
the form D + A (I) for some effective divisor D, we have natural linear maps 

g m : I m H° (f, O r (mH - A(/))) , 

for all m > 0. 


Remark 3.1. Consider the exact sequence 

0 —y I Or —y 7r* (IOjI) —y J~ —y 0, 

where I is the ideal sheaf associated to I, and T is the cokernel. Taking global 
sections, we get, for m 3> 0, the exact sequence 

0 —y H 0 (r,IO r (m)) —y H°(T,IO T (mH )) —y H° (T,P) -y 0. 

^The term point will always refer to a closed point. 
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Indeed, T has finite support and, since the normalization map 7r is finite, we have 
H°(TjO r (mH)) 9* H°(T, 7r*(I0p)(m)). Since IO r (mH) = 0 T (mH - A(/)) and, 
for m S> 0, H°(V, IOr(m)) = / m //(r) m , we get, for m 3> 0, the exact sequence 

0 -A 7 m /7(r) m W 77° (r, O r (m77 - A(7))) -a 77° (r, T 7 ) -A 0. 

In particular, for m > 0, 

ker(£ m ) = 7(r) m . 

Definition 3.2. With notation and assumptions as above, the ideal 7 is called an 
adjoint ideal of T if the maps 

Qm : I m -A 77° (r, Op ( mH - A(7))) 

are surjective for m large enough. 


As already remarked in the introduction, the existence of adjoint ideals is classical. 
Locally, adjoint ideals are characterized by the following criterion: 

Theorem 3.3. The ideal I is an adjoint ideal of T iff Ip = IpOr,p f or Cl P £ 

Sing(r). 


Proof. Using the notation of Remark |3.1[ we have, for m 3> 0, 

dim,, coker g m = h° (r, P) = ^ l(IpO T ,p/Ip). 

PeSing(r) 


Hence, Q m is surjective iff IpOr,p = Ip for all P 6 Sing(r). 

Corollary 3.4. If I is an adjoint ideal ofT and P £ Sing(r), then Ip C Or,p- 
Proof. Suppose Ip = Or,p- Then Ip C IpOr^p, a contradiction to Theorem 
Corollary 3.5. The support of Proj (S/1) contains Sing(r). 

Proof. Follows immediately from Corollary |3.4[ 


3.3 


□ 


□ 


□ 


Theorem 3.6. There is a unique largest homogeneous ideal 0 C S which satisfies 

0P = Co r P for all P £ Sing(r). 

The ideal 0 is an adjoint ideal of T containing all other adjoint ideals of T. In 
particular, 0 is saturated and Proj(5/0) is supported on Sing(r). 

Proof. For the conductor ideal sheaf C = Ann£) r (7r*Ctp/C>r) on T, we have Cp = 
C-o r p f° r P £ r. If j : T —» P/ is the inclusion, then the graded 5-module 
0 = (B ne zH 0 (F r k , j*C(n)) associated to j*C is the unique largest homogeneous ideal 
with 0p = Co r P for all P £ Sing(r). By Theorem 


3.3 


and the properties of the 
conductor, 0 is an adjoint ideal. Moreover, if I is any other adjoint ideal, then 
lpC0p for all P £ T, hence I C 0. □ 


Definition 3.7. With notation as in Theorem 3.6, the ideal 0 is called the Goren- 
stein adjoint ideal of T. We also write 0(r) = 0. 

For repeated subsequent use, we introduce the following notation: 


Notation 3.8. Let T C Pt be a curve as above with Gorenstein adjoint ideal 0. Let 
C be the affine part of T with respect to the chart 

A l -A P[., (Ad,..., Ad) ha- (1 : Xp : • • • : X r ) 

let 1(C) C k [X u ..., X r \ be the ideal of C, let k[C\ = k[x \,..., x r ] = k [Ad,..., X r \ /1(C) 
be its coordinate ring, and let Sing(C') be its set of singular points. 
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Proposition 3.9. Assume T C PI is a curve as in Notation 3.8, with affine part 
C. Let (5 be the ideal of k[C\ obtained by dehomogenizing © with respect to Xo and 
mapping the result to k[C]. Then 


® — Cfc[C] ■ 

If r has no singularities at in/frora 2 and C^c\ = {Si( x i> ••• ,x r ) | i)k\c\ w ' 1 ^ 1 polyno¬ 
mials gi G k[X i, ...,X r \, then © is the homogenization of 


<^(Xi,..,X r )|i) fc[Xli ... iXr] +/((7) 

with respect to Xq. 


Proof. The first statement is obtained by localizing at the points of C: 

©p = Co c ,p = ( C k[C\)p for eac]l P £ C. 

Here, the first equality is clear from the definition of © (see Theorem 3.6). The 
second equality holds since forming the conductor commutes with localization since 
k[C] is normalization-finite (see |Zariski and Samuel 1975, Ch. V, § 5]). 

The second statement of the proposition follows from the first one since there are 
no singularities at infinity, © is saturated, and the support of © is contained in C. 

□ 


We take a moment to specialize to plane curves. 


Remark 3.10. Assume T is a plane curve. Then, by Max Noether’s Fundamentalsatz, 
the maps g m : © m —> H° (T, Op {m.H — A(©))) are surjective for all m. Referring 
to each homogeneous polynomial in © not contained in I(T) as an adjoint curve to 
T, this means that residual to A(©), the adjoint curves of any degree m cut out the 
complete linear series A m = \mH — A(©)|. See [van der Waerden 1939, § 49]. 

Theorem 3.11. Assume T is a plane curve of degree n. Then, residual to A(©) ; 
the elements of © n _3 cut out the complete canonical linear series. Equivalently, 


deg A(0) = 2«5(T). 


Proof. See [Gorenstein 1952 


Thm. 9], 


(3.1) 

□ 


Recall that the dimension of the canonical linear series is dimM n _3 = p(T) — 1. 


Remark 3.12. Assume T is a plane curve of degree n. If p(T) = 0, that is, T is 
rational, then dim„4 n _2 = deg A n -2 = n — 2. In this case, the image of T under 
A n - 2 is a rational normal curve V n -2 © P£~ 2 of degree n — 2. Via the birational 
morphism T n _2 —> T, the problem of parametrizing T is reduced to parametrizing 
the smooth curve T n - 2 - For the latter, we may successively decrease the degree 
of the rational normal curve by 2 via the anti-canonical linear series. This yields 
an isomorphism from V n -2 either to P 1 or to a plane conic, depending on whether 
n is odd or even. If T is defined by an equation over a subfield L C k, then all 
computations considered so far take place over the coefficient field L. Parametrizing 
the conic, however, may require a quadratic field extension, depending on whether 
the conic contains an L-rational point or not. See Bohm 1999 and Bohm et al. 
2015c| for details. 


By generalizing the formula in Theorem 3.11 
of adjoint ideals 


we now derive a characterization 


which is also valid in the non-planar case. We use the following 
notation: If I C S is a homogeneous ideal, write deg I = deg Proj (S/I). That is, 
deg / is (dim / — 1)! times the leading coefficient of the Hilbert polynomial of S/I. 


9 

If k is infinite, this assumption can always be achieved by a projective automorphism defined 
over k. Otherwise, we may have to replace k by an extension field of k. 
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Lemma 3.13. Let I C S be a saturated homogeneous ideal with I(T) C I. Then 

deg A (I) < deg I + S(T), 
and I is an adjoint ideal of T iff 

deg A(7) = deg I + <5(T). 


Proof. Let Pr(i) = (degT) • t — p a (T) + 1 be the Hilbert polynomial of fc[r]. Denote 
by Jr the image of / in fc[r]. Then, for m 0, 


deg I = dim k (S m /I m ) = dim fc (7[r] m /(/r) m ) = Pr(m) - dim fc (7 r ) m . 

h° (r, O r (mH - A(7))) = dirn fc (7 r ) m + h° (T,T) > 


3.1 


Moreover, by Remark 
dimfc(/r)m for m ^$> 0. Hence, by Riemann-Roch, for m > 0, we have 


(degT) • m — deg A (I) = deg \mH — A(/)| = dim \mH — A(/)| +p(T) 

> dinifc (I r )m - 1 +p{T) 

= Pr(m) — deg I — 1 +p(r) 

= (degT) • m — <5(r) — deg I. 


Here, we use that |mi7 — A(/)| is nonspecial for large m by reason of its degree. 
Equality holds iff g m is surjective. □ 


Remark 3.14. In the case where T is a plane curve and Z = (25 is its Gorenstein 
adjoint ideal, Lemma 3.13 shows that Equation ( |3.1[ ) may be rewritten as 

deg © = <5(r). (3.2) 


Note that (3.1) and (3.2) may not hold in the non-planar case: 


Example 3.15 ( de Jong and Phster 2000, Example 5.2.5]). Let T C 
image of the parametrization 


be the 


ji 
■ € 


(s : t) i-A (s 5 : fV : f 4 s : t 5 ). 


Then T has exactly one singularity at (1 : 0 : 0 : 0). Furthermore, p(T) = 0 and 
p a (r) = 2, hence J(r) = 2. However, C5 = {X\, X2, Xf) C C[Xo,..., A 3 ], hence 
deg © = 1 . 


Remark 3.16. If T C Pj] is any curve as in Notation 3.8, with affine part C and no 
singularities at infinity, then it follows from Proposition |3.9| that 


deg© = dim fc (k[C\/C k[c] ) = EpeSmgRO dim k (O c ,p/C OC}P ). 


Lemma 3.17. If char k = 0, then dinifc(C>r,p/Co r p) ^ 6 p(r) for any point P & T. 


Proof. Since normalization commutes with base change, this follows from the case 
k = C proved in |Greuel 1982, 2.4], □ 

Now recall that a point P S Sing(r) is called a Gorenstein singularity if 

dinu, (Or,p/Cc>r,p) = ^pCO- 


Example 3.18. Plane curve singularities are Gorenstein (see, for example, jde Jong] 


and Phster 2000 Corollary 5.2.9]). 


Corollary 3.19. We have: 

(1) If char k = 0, then deg© < <5(r). 

(2) If r has only Gorenstein singularities, then 

deg© = J(r) and deg A(©) = 25(r). 


Proof. This is clear from the discussion above. 

□ 
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In the case of arbitrary singularities, we will make use of the equality 

deg © = deg A(0) — <5(T) 

to compute deg © without actually knowing ©, and apply this in the final verification 
step of our modularized adjoint ideal algorithm. To this end, if k = C and T is defined 
over the rationals, we will present a modular approach to computing deg A(©), and 
we will use standard techniques to compute <5(T). In fact, for the latter, first note 
that the delta invariant of T differs from that of a plane model of T by the quantity 
p a (T) — ( deg ^ -1 )- The delta invariant of a plane curve, in turn, can be computed 
locally at the singular points, either from the semigroups of values of the analytic 
branches of the singularity (see [de Jong and Pfister 2000| , Greuel et al. 2007] ), or 
from a formula relating the local delta invariant to the Milnor number (see Remark 
Section [ 7 ] below). 


7.3 


m 


Remark 3.20. Note that computing deg A(©) also mean s to compute the dimension 
Given T 


dim k{k[C]/C k [cj): Given I C 
singularities at infinity, we have 


as in Notation 


with affine part C and no 


deg A(©) = <5(T) + deg© 


= dim*. (k[C]/k[C]) + dim*. (k[C}/C k[c] ) 

= dim*. (k[C\/C k[c] ). 

We are now ready to address the computation of the Gorenstein adjoint ideal. 
Using Pro posit ion |3.9[ one way of finding © is to apply the global algorithm presented 


in Section 4.2 below, starting from the normalization k[C\. The normalization, in 


turn, can be found by combining the minimal local contributions to k[C] at the 
singular points via Proposition 2.11 As it turns out, however, it is more efficient to 


directly compute local Gorenstein adjoint ideals at the singular points, and get © 
as their intersection. This will be discussed in Sections [5]and[6j 


Remark 3.21. In applications, T C is often defined over a perfect subfield k' C k 
(for example, k! = Q and k = C). In such a situation, by base change, <5(T) = 
6 (T(k')). Moreover, since the algorithms in Sections [5] and [6] rely on Grobner bases, 
and Buchberger’s algorithm for computing Grobner bases does not leave the ground 
field, ©(T) = &(T(k'))K[X 0 ,..., X n ], and generators can be found by computations 
over k'. 


4. Global approaches 


4.1. Computing the conductor via the trace matrix. We will require some 
facts from classical ideal theory (see |Zariski and Samuel 1975, Ch. V] for details 
and proofs): Let R be an integral domain, and let K = Q (R) be its quotient field. A 
fractionary ideal of R is an R-submodule b of I\ admitting a common denominator: 
there is an element 0/d£l2 such that d b C R. 


Example 4.1. The extensions Aj computed by the normalization algorithms from 
Section [2] are fractionary ideals of the given affine domain A. 

If b, b' are two fractionary ideals of R, with b 1 nonzero, then b : b' = {z 6 K \ 
z b’ C b} is a fractionary ideal of R as well. A fractionary ideal b of R is invertible 
if there is a fractionary ideal b ’ of R such that b • b' = R. In this case, b 1 is uniquely 
determined and equal to R : b. 

Suppose in addition that R is normal. Let K' be a finite separable extension of 
K, and let R' be an integral extension of R such that K' = Q {R')- Moreover, let 

T vk'/k ■ K ' K, zia ^2 d( z ), 

gSGal(A''/A') 
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be the trace map. Then the complementary module 

£r'/r ■= {z E R' \ Tr K ,/ K ( zR') C R} 

of R! with respect to R is a fractionary ideal of R' containing R'. Hence, the different 


'Sr'/r — R' '■ £r'/r — { z ^ R' 


z£ 


R'/R 


c R!) 


= {z£ K' I zx E R' for all x E R' with Tr k'/k (xR') C R} 
of R' over R is a nonzero ideal of R!. 

Now, keeping our assumptions, we focus on the case where R is a Dedekind 
domain, and where R' is the integral closure of R in K'. Then R' is a Dedekind 
domain as well, which implies that every nonzero fractionary ideal of R' is invertible. 
On the other hand, by the primitive element theorem, there is an element y E R' 
with K' = K(y). Denote by f(Y) E K[Y] the minimal polynomial of y over K. 
Then, as shown in |Zariski and Samuel 1975, Ch. V], 

f'{y)R' = CRI/R[y]T>Rf /R, 

hence 

Cr'/R[v\ = f\y)^R'/R- l 4 - 1 ) 

We now fix the following setup: 


Notation 4.2. Let k be a perfect field. Let T C be a plane curve of degree n defined 
by an irreducible polynomial F E k[X, Y , Z\. Suppose that T has no singularities at 
infinity with respect to the affine chart 

A (X,Y)^(l:X:Y), 

and that the equation / E k[X, Y] of the affine part C of T is monic in Y. 


Write k[C\ = k[x,y\ = k[X,Y]/(f(X,Y)) for the affine coordinate ring of C and 
k(C) = k(x,y) = k(X)[Y\/(f(X,Y)) 

for its function field. Then x is a separating transcendence basis of k{C) over k, 
and y is integral over k[x], with integral equation f(x,y) = 0. In particular, k[C] is 
integral over k[x\, which implies that k[C\ coincides with the integral closure k[x\ of 
k[x] in k(C). Furthermore, k[C] is a free A;[x]-module of rank 

n ■= deg y (f) = [k(C ) : k(x)}. 


Definition 4.3. An integral basis for k[C\ is a set bo,, b n -± of free generators for 
k[C\ over k[x\. 


k[C] = k[x]bo © • ■ ■ © k[x]b n -\. 


Remark 4.4. Since k(C) = k(x,y) = k(X)[Y]/(f), any element a E k(C ) can be 
represented as a polynomial in k(X)[Y] of degree less than n = deg/. Hence, 
one can associate to a a well-defined degree deg y (a) in y and a smallest common 
denominator in k[x] of the coefficients of a. In particular, k[C] has an integral basis 
(bi) in triangular form, that is, with deg y (bi) = i, for i = 0, ...,n — 1. If not stated 
otherwise, all integral bases will be of this form. In principle, such a basis can be 
found by applying one of the normalization algorithms discussed earlier. However, in 
the characteristic zero case, methods relying on Puiseux series techniques are much 
efficient (see |Bohm et al. 2015a| 


more 


and 


van 


Hoeij 1994]). 


Example 4.5. An integral basis for the curve from Example 2.10 is given below: 

y(y- i) 2/(2/ — i ) 2 y 2 (y- i ) 2 


L 


y, 


3.9 


Using Proposition 
and K' = k(C), we get Algorithm [T] 


and Equation (4.1), with R = k[x \, R! = k[C\, K = k(x), 
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Algorithm 1 Gorenstein adjoint ideal via linear algebra (see 


Mnuk 1997 ) 


Input: A plane curve T with affine part C as in Notation |4.2| 
Output: The Gorenstein adjoint ideal © of T. 

1: Compute an integral basis {bi ) i=0 n _ l for k[C]. 

2: Compute the (symmetric and invertible) trace matrix 


T = OIAioaw m).._ 0 .6 *(*)"*”. 


3: Compute a decomposition L ■ R = P ■ T, where L is left triangular matrix 
with diagonal entries equal to one, R is a right triangular matrix, and P is a 
permutation matrix. 

4: For j = 0 ,n — 1, use forward and backward substitution to compute 


n —1 

r h = 

i=0 


Sij bi , 


where () = T 1 . The rjj are /c[x]-module generators for ^ucyk \ x ]• By (4-1), 
Ck{C] = (|y ( x ^v)Vj I j = 0,...,n — 1). 

5: Let C be the ideal of k[X,Y] generated by representatives of minimal y-degree 
of the f£ ( x,y)r]j, j = 0, 1. 

6: return the homogenization of C with respect to Xq. 


Remark 4.6. To compute an integral basis via Puiseux series in the characteristic 
zero case, we temporarily may have to pass to an algebraic extension field of k. 


Example 4.7. The curve TCP?, from Example 2.10 with affine equation 


X 5 - Y 2 (1 - Y ) 3 = 0 


has a singularity of type A 4 at (0,0) and a 3-fold point of type Eg at (0,1). From 
the integral basis 


6 0 = 1, h=y, b 2 = v(v x , b 3 = , and b A = 

given in Example |4.5[ we compute the trace matrix 


T = 


( 5 

3 
0 
0 


3 

3 

0 


0 0 0 \ 

0 0 —5x 2 

0 — 5x 2 —3x 


0 -5x 2 —3x 


0 

0 J 


\ 0 -5x z —3x 0 

which yields by forward and backward substitution 

Cc[C\ = (■ x3 , x 2 (y- 1), xy (x - 1) , y{y- l) 2 ^ 


C[C] 


Homogenization (and primary decomposition) gives 

© = (y, x 2 ) n (|a 3 , x(y-z),(y- z) 2 ^ . 

Writing © as the intersection of simpler ideals corresponding to the singularities 
motivates the local to global approach discussed in Sections [5] and [6] below, where 
© will be found as the intersection of local Gorenstein ideals. 
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4.2. Computing the adjoint ideal via ideal quotients. The algorithm pre¬ 
sented in what follows relies on normalization and ideal quotients. It is not limited 
to plane curves. 

Proposition 4.8. Let T C Pf be a curve with affine part C as in Notation 
Write k[C] = \U, where U C k[C] is an ideal and d E U is nonzero. Then the 
conductor is 

Ck[C] = {d) k [c} '■ u ■ 

Proof. By definition, 

C k[c] = {sek[C]\s-k\C]ck[C}} 

= {s£ k[C] | s ■ g £ (d) for all g E 17} 

= (d)k[C] : U - □ 

Using once more Proposition |3.9[ we get Algorithm [2] 

Algorithm 2 Gorenstein adjoint ideal via ideal quotients 

Input: A curve T C Pf with affine part C as in Notation 
at infinity. 

Output: The Gorenstein adjoint ideal 0 of T. 
l: Normalization: Compute polynomials d, ao,... ,a s E k[X i,..., X r \ such that the 
fractions generate k[C] as a k[C] -module. 

2: Compute the ideal quotient 

C = ((d) + 1(C)) : ((o 0 ,...,o s ) + /(C)) C k[X 1 ,...,X r }. 

3: return the homogenization of C with respect to Xq. 


3.8 and no singularities 


3.8. 


Example 4.9. In Example |2.10[ 

a 0 = X 3 , a v = X 2 Y(Y - 1), a 2 = XY(Y- l) 2 , a 3 = Y 2 (Y - l) 2 , 
and d = X 3 . Hence, 

(d, f) : (ao, o 3 , /) = (X 3 , X 2 (Y - 1), XY (Y - 1), Y (Y - l) 2 ) . 

5. A Local to global Approach 

In this section, motivated by the local to global approach for normalization, we 
introduce local Gorenstein adjoint ideals of a given curve and show how to find 
the Gorenstein adjoint ideal 0 as their intersection. Together with the algorithm 
presented in the next section, where we will show how to compute the local ideals, 
this yields a local to global approach for finding 0. As we will see in Section [9j 
this approach is per se faster than the algorithms discussed so far. In addition, it is 
well-suited for parallel computations. 

We fix the following setup: 

Notation 5.1. Let T C P£ be an integral non-degenerate projective curve, and let S 
be the homogeneous coordinate ring of PL 

Definition 5.2. Let W C Sing(T) be a set of singular points of T. The local 
Gorenstein adjoint ideal of T at W is defined to be the largest homogeneous ideal 
0(W) C S which satisfies 

<S(W) P = C 0v p for all P E W. 

For a single point P E Sing(T), we write 0(P) := 0({P}). 


(5.1) 
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Remark 5.3. Since ©(VP) is the largest homogeneous ideal satisfying (5.1), it is 
saturated and Proj(5/0(W)) is supported on W. 


Proposition 5.4. Let W C Sing(r). Then 


&(w) = n P£W &(p). 

Proof. This is immediate from the definition: If 05' := (~) PgW 05(P), then Proj(S'/© / ) 
and Proj(«Sy©(VP)) have the same support W, and 


<5' Q = 05(Q) Q = Co riQ = <5(W)Q 
for all Q E W, hence ©(VP) = ©'. 


□ 


Proposition |5.4| yields Algorithm [3j 


Algorithm 3 Gorenstein adjoint ideal, local to global 


as in Notation 5.1 


Input: A curve T C 
Output: The Gorenstein adjoint ideal 05 of T. 
1: Compute Sing(T) = {Pi, ...,P S }. 


Apply Algorithm [4] in Section [b] below to compute © (Pf) for all 
return f)i=i® ( P i)- 


i. 


Remark 5.5. It is clear from Proposition 5.4 that we may choose any partition 
Sing(T) = Uf=i Wi °f Sing(T) and have 

© = 

This is useful in that for some subsets Wi, specialized approaches or a priori knowl¬ 
edge may ease the computation of ©(IT*). In Section [T| we will present some ideas 
in this direction for plane curves. 


6. Computing local adjoint ideals 


In this section, we modify Algorithm [2] so that it computes the local Gorenstein 
adjoint ideal at a point P from a minimal local cont ribution at P via ideal quotients. 

with affine partj^] C and a point 


We consider a curve T C PI as in Notation 


3.8 


P E Sing(C'). Let be the minimal local contribution to k[C\ at P; so U C k[C\ 
is an ideal and d E U is nonzero. 

Proposition 6.1. With notation as above, and given Q E C, we have 

Coc o tf Q = P, 


((d) k [c] 

Proof. By the minimality assumption, we have 


■■U)q = { 


Oc’q if Q + P- 


U = 


_ O c ,Q if Q = P, 


1 

d L J Q \ O c ,q if Q + P- 

The claim follows since localization commutes with forming the conductor: 


k[C] 


:U) = (Cm 


q ~ \^(- d u)/m) Q - c (\u) Q / k [c] Q 


□ 


^To cover all singular points of T, we may have to choose affine charts other than that considered 


in Notation 


3.8 
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Now, we argue as in the proof of Proposition |3.9 
Remark 


5.3 


From Proposition |6.1| and 


it follows that (d) k i C i : U coincides with the ideal obtained by deho- 
mogenizing 0(P) with respect to Xq and mapping the result to k[C\. Hence, since 
0(P) is saturated, Algorithm [4] below indeed computes 0(H). 


Algorithm 4 Local Gorenstein adjoint ideal from local contribution 


Input: A curve T C P7 with affine part C as in Notation 
Sing(C). 

Output: The local Gorenstein adjoint ideal 0(P) of C. 
l: Compute polynomials d,ao,...,a s E k[X i,...,A r ] such that the fractions 
< 1 < 1 (t i 'x r ) g enera t e the minimal local contribution to k [C] at P as a k [ C ]- 
module. 

2: Compute the ideal quotient 


3.8 and a point P E 


C = ((d) + I (C)) : ((a 0 ,..., a s ) + I (C)) C k [X u ...,X r ] . 


3: return the homogenization of C with respect to Xq. 


Example 6.2. We compute the local Gorenstein adjoint ideals for the curve given 
in Example |2.10 with affine equation 

X 5 - Y 2 (1 - Y ) 3 = 0. 

For the AG-singularity Pi, we found 

di = x 2 and U x = < x 2 , y(y - 1) 3 > C[C] , 

so that 

0(Pi) = (A 2 , Y) . 

For the Eg singularity Po, we observed that 

d 2 = x 3 and U 2 = (x 3 , x 2 y 2 (y - 1), y 2 (y - l) 2 ^) 

leading to 

0(P 2 ) = (A 3 , A (Y - Z ), (Y - Zf) . 


€[C] 


Note that 0(Pi) and 0(P 2 ) are the ideals already obtained in Example 4.7 


7. Improvements to the local strategy for plane curves 
In this section, we focus on the case of a plane curve F with affine part C = V (/) 


and Sing(r) = Sing(C) as in Notation 4.2 For simplicity of the presentation, we 


suppose throughout the section that our ground field k = C. 

As explained in Section [5j the Gorenstein adjoint ideal 0 can be computed as 
the intersection of local Gorenstein ideals via a partition of Sing(C). To begin with, 
consider the following partition: 


Sing(C) = W 2 U W 3 U • • • U W r U W', 


(7.1) 


where, for all i, W t denotes the locus of ordinary Pfold points (ordinary multiple 
points of multiplicity i), and where W' collects the remaining singularities of C. In 
particular, W 2 is the set of nodes of C. 

Lemma 7.1. Let P E Sing(C'), and letmp C k[X,Y] be the corresponding maximal 
ideal. If P is an ordinary i-fold point of C, then 


&(P) = m'p 


i— 1 
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Proof. Since C is a plane curve and P is an ordinary z-fold point o f C, t he cond ucto r 


Cq cp = m l C p, where m c,p is the maximal ideal of Oc,p (see 


Maths 1970 


and Valabrega 1979]). The result follows from the very definition of ©(P). 


Greco 


□ 

Applying the lemma to the partition S we get the intersection of ideals 

© = / (w 2 ) n i (w 3 ) 2 n • • • n / ( w r ) r n < 3 {w'). (7.2) 

Hence, in the case where T is known to have ordinary multiple points as singularities 
only (that is, W' = 0), we can compute <25 in a very efficient way by using Algorithm 
[5] below (see Bohm 1999]). 


Algorithm 5 Gorenstein adjoint ideal, ordinary multiple points only 


Input: A plane curve T of degree n with defining polynomial F as in Notation 4.2 
with only ordinary multiple points. 

Output: The Gorenstein adjoint ideal <25 of T. 


T •— / dP dP. dF\ 
\dX- dYi dZ / 

i := 1 

while (J t : (A, Y, Z) c 
i := i + 1 

j I d J+l+m F 
Jl ■ \axWY l dz m 

B := (X, Y. Z) n_i 

while i > 0 do 


(the ideal defining Sing(T)) 

3 ) / (1) do 

j + l + m = i, j, l, m £ No 


8 

9 

10 

11 


h '■= (Ji -1 : B°°) (the ideal of the z-fold points of T) 

B := ((B Pi l\~ l ) : (X, Y, Z)°°) 
i := i — 1 

return B 


In the general case, Equation (7.2) allows us to reduce the computation of & to 
the less involved task of computing <3{W') as soon as we detect the ordinary z-fold 
points. To begin with treating these, here is how to find the nodes: 

Remark 7.2. We know how to find all singularities: Sing(C') is given by the ideal 


J = 


dX' dY 


By the Morse lemma (see Milnor 1968]), a point P E Sing(C) is a node iff the Hes¬ 
sian matrix Hess(/) formed by the second partial derivatives of / is non-degenerate 
at P. That is, P is a node iff 

I(P) + (det(Hess(/))) = k[X,Y]. 

This gives us a fast way of computing W 2 . 


Carrying our efforts one step further, we discuss the local analysis of the singu¬ 
larities via invariants. This yields an efficient method not only for finding the delta 
invariant, but also for detecting the ordinary z-fold points, for each z: 


Remark 7.3. Let P E Sing(C'). After a translation, we may assume that P = (0,0) 
is the origin. Write mp for the multiplicity and 


= dim*. 



d.f df\\ 
dX' dY/) 


for the Milnor number of C at P. Then mp = deg/i p , where hp is the lowest 
degree homogeneous summand of the Taylor expansion of / at P. Recall that hp 
can be computed via standard bases (see |Greuel and Pfister 2008 ). Furthermore, 
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if the Newton polygon of / is non-degenerate (otherwise, successively blow up), the 
number of branches of / at P can be computed as 


£”> d 


r P = 


V' 


O+i) 


x 


- V, 


0) y 0+i) 


x 


Y 


~ \\ 


0 ) 


Y 


where RW,I/( s ) are the (ordered) vertices of the Newton polygon (and X and Y 
refer to the respective coordinates). This is immediate from |Brieskorn and Knorrer 


1986, Section 8.4, Lemma 3]. The delta invariant of C at P is then obtained as 


5p — -(gp + r P — 1 ) 


(see, for example, Greuel et al. 2007, Chapter 1, Proposition 3.34]). Furthermore, 
P is an ordinary z-fold point iff hp is square-free and mp = i. Equivalently, 


(■ m P ,r Pl 5p) = [i, i 


See Greuel et al. 2007, Chapter 1, Proposition 3.33]. 


The local analysis of the singularities may be used to further refine our partition 
of Sing(C). For example, singularities of type ADE can be identified as follows: 


Remark 7.4. With notation as in Remark 7.3, the point P = (0,0) £ Sing(C) is 

(1) of type A n , n > 2, iff hp = l\, with l\ £ k[X,Y] linear, and /ip = n, 

(2) of type D n , n > 4, iff hp = I 1 I 2 I 3 ° r hp = Zf Z 2 , with pairwise different linear 
polynomials lj £ k[X,Y], and gp = n, and 

(3) of type E n , n = 6, 7, 8 , iff hp = if, with l\ £ k[X, Y] linear, and gp = n. 
Here, in (2), hp splits into three different linear factors iff P is of type ZL 4 . See, 

for example, [Greuel et al. 2007], Chapter 1 , Theorems 2.48, 2.51, 2.54], 

To describe the local Gorenstein adjoint ideal at a singularity of type A. D , or E, 
we use the following notation: 


Notation 7.5. For any element g £ k[[X, Y"]], let gj = taylor (g,j) £ k[X,Y] be the 
Taylor expansion of g at P = (0,0) modulo 0(j + l)j^] 


If C has a singularity of type A n at P = (0,0), we may write / in the form 
/ = T 2 + W n+1 , where T, W £ k[[X, T]] is a regular system of pa ram eters. Let 
s = (the meaning of s will become clear in the proof of Lemma 

compute the Taylor expansion T s _i £ k[X,Y] as follows. If n and t 
to 1, set Tq = 0. Otherwise, inductively solve / for T: 


7.6). We may 


thus s is equal 
Start by choosing a linear 


form Ti £ k[X, Y] such that taylor(/, 2) = T 2 . Supposing that 1 < j < s — 1 and 
Tj = T + 0(j + 1) has already been computed, write 

taylor(/ - T 2 , j + 2) = 2 T x ■ m, 

with m £ k[X , Y] homogeneous of degree j + 1, and set Tj + \ = Tj + m. 


Lemma 7.6. Let C have a singularity of type A n , n > 1, at P = (0,0). Set 
s = and let T s _\ be defined as above. Then C5(-P) is the homogenization of 

(X s , T a _!, Y s )ck[X,Y] 

with respect to Z. 


Proof. The case n = 1 is clear, so we may suppose n > 2. If & = (X s , T s _ 1 , W) C 
k[X, Y], then &q = Oc,q for all Q £ C \ {T}, so it suffices to show that & p = Cp, 
where B = Oc,p • For this, we pass to the completion 

B = k[[x,y]} = k[[X,Y}}/(f(X,Y)), 

4 The notation 0(m ) stands for terms of degree > m. 
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and consider the isomorphism 

A = k[[t, w}] = k[[T, W]]/ (T 2 + W n+1 > T(x, y),w ha W(x, y ). 

An analysis of the normalization algorithm applied to A shows that 


^ = E*[W]-» i + E fc iwi • Ti 


i=0 


i=n—s +1 


Bohm et al. 2014 


Sect. 4]). 


and that it takes s = steps to reach A (see 

Hence, 

C A = (t, w s ) A , so that Cg = (T(x, y ), W (. x , y) s )g- 
Working in k[[X, Y - ]], we write 

T = aX + bY and W = cX + dY, 

with a, b,c,d& k[[X, y]] and such that ad — be is a unit in k[[X, y]]. Since (X, Y) = 
(T, W), it follows that (X,Y) S = {T,W) S C ( T,W S ). Since (X,Y) = (X,T) or 
(X,Y) = (T, Y), we have W s € (X,Y) S C (X S ,T,Y S ). We conclude that 

(X S ,T,Y S ) = (T, W s ). 

If s > 1, then (X , Y) = (X, T s _i) or (X, Y) = (T s _i, Y), hence, for any s, we have 
(X, Y) s C (A s , T s _i, Y s ). We conclude that 

(X S ,T S _ 1 ,Y S ) = (X S ,T,Y S ). 

Now recall that B is an excellent ring, which implies that B = B (see, for example, 
[Bohm et ah 2014, Sect. 1]). It follows that 

Cg = Homg (^B, B'j = Horns (B, B) ®b B = Cb ®b B. (7.3) 

Since completion is faithfully flat in the case considered here, we conclude that 

C B = {x s ,T s _i(x,y),y s ) B . 

□ 


Remark 7.7. In particular, if P is a cusp, then &(P) 
and cusps may be treated simultaneously. 


(X, Y). So, in (7.2), nodes 


If C has a singularity of type D n at P = (0,0), we may write / in the form 
f = W ■ (T 2 + W n ~ 2 ), where T,W € k[[X, Y"]] is a regular system of parameters. 
Let s = |_|J. We may compute the Taylor expansion T s _2 G k[X,Y] as follows. 
If n = 4, set To = 0. If n > 5, choose linear forms T\ , W\ G k[X,Y\ such that 
taylor(/, 3) = Tf ■ W\. For j < s — 2, determine Wj = W + 0(j + 1) as the Puiseux 
expansion up to order j of / corresponding to W\. Supposing that 1 < j < s — 2 
and Tj = T + 0{j + 1) has already been computed, write 

taylor(/ - T 'j ■ W j+ 1 , j + 3) = 2Z X -Wi-m, 

with m G k[X, Y] homogeneous of degree j + 1, and set Tj +1 = Tj + m. 

Lemma 7.8. Let C have a singularity of type D n , n > 4, at P = (0, 0). Set s = |_^J, 
and let T s _2 be defined as above. Then &(P) is the homogenization of 

(X, Y) ■ (X s - 1 , T s _ 2 , y*- 1 ) C k[X,Y] 


with respect to Z. 
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Proof. We have an isomorphism 

A-$ B, q i-a T(x,y), w i-)- W(x,y), 


where B = Oc,p and 

A = k[[t, w]] = k[[T, W]}/( W • (T 2 + W n ~ 2 )) . 
This time, the normalization is 


n—2—s 


n —3 


a= E MW]•w i + E mw]• t+*[[*]] 


w 


n—2 


2=0 


i=n—l—s 


and it takes s = |_§J steps to reach A (see again 


t ' t 2 ’ 

Sect. 4]). Hence, 


Bohm et al. 2014 


Write 


Ca = (t 2 , tw, w 3 ) . 

T = aX + bY W = cX + dY 


with a, b,c,d£ k[[X, Y]] and such that ad — be is a unit in k[[X, Y]]. Since ( X , Y) = 
(T,W), we have (XT, YT) = (T 2 ,TW) and ( X,Y) S = (T,W) S C {T 2 ,TW,W s ), 
hence 

(X,Y) ■ (X s - 1 , T,^” 1 ) C (T 2 ,TW,W s ). 

For the other inclusion, observe that (X,Y) = (X, T) or (X,Y) = (T,Y ), so it 
follows that (I,h ) s_1 C (X s ” 1 , T,Y s ~ l ), hence 

W s € (X,Y) S C (X,Y) ■ (X s - 1 , T,^” 1 ) . 

If s > 2, then (A", Y) = (X, T s _ 2 } or (A, Y) = (T s _ 2 , Y), hence, for any s, we have 
(A, Y ) s-1 C (A s - 1 ,T s _ 2 ,Y' s " 1 >. We conclude that 

(A*” 1 , T s _ 2 , Y 3 - 1 ) = (X s - 1 , T, Y 3 - 1 ) . 


To summarize, 

(T 2 ,TW,W s ) = (A, Y) ■ (X s - 1 , T^ 3 - 1 ) = (A, Y) • (X s ' 1 , T s _ 2 , Y 5 " 1 ) , 

hence 

C# = <ar,2/> • (.T s_ 1 ,T s _ 2 (x,y),y s_1 ) C B. 

Then the claim follows as before. □ 


Lemma 7.9. Let C have a singularit y of type E n , n = 6,7,8, at P = (0,0). Set 
s = ; and let l\ be as in Remark l.f. Then C5 (P) is the homogenization of 


(A, Y) ■ (A s_1 , h, Y s_1 ) C k[X, Y] 


with respect to Z. 


Proof. Depending on n E {6, 7, 8}, we have an isomorphism 

A-* B, q i-a T(x,y), w i-)- W(x,y), 

where B = Oc,p and 

A = k[[t,w]] = k[[T,W}]/ (T 3 + W A ), 

A = MM]] = k[[T, W]]/ (T (T 2 + W 3 )), 
A = k[[t,w]} = k[[T,W]}/(T 3 + W 5 ), 
respectively. In each case, by (Bohm et al. 2014, Sect. 4], 

A = k[[w]}-l + k[[w]}-^ + k[[w}}~, 

which implies that 

Ca = (t 2 ,tw, w 3 ') . 
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The same argument as in the proof of Lemma |7.8| shows that 

Cb = ( X ’V) ' (z s_1 > T s-2(ahy),2/ s_1 ) C B, 
and the claim follows as before. Note that T s _2 = 0 if s = 2, and T s _2 
s = 3. 


h if 

□ 


In principle, we could pursue a similar strategy for all singularities classified by 
Arnold in Arnold et al. 19951. However, in Bohm et al. 2015a], we give an 
algorithm which, for plane curves in characteristic zero, allows us to compute the 
local contributions to the normalization for a broad class of singularities in a direct 
way. Combining the approach of Section [6] with this algorithm or with modular 
techniques and normalization as described in Section [8] below, we already get a very 
efficient algorithm for computing ®. 


Remark 7.10. For the local analysis of the singularities, we temporarily may have 
to leave k. 


8. Parallel computation and modular techniques 


Algorithm [3] is parallel in nature since the computations of the local adjoint ideals 
do not depend on each other. In this section, in the case where the given curve is 
defined over Q, we describe a modular way of parallelizing Algorithm [3] even further. 
One possible approach is to replace the computations of the Grobner bases involved, 
the computation of the (minimal) associated primes in the singular locus, and the 
computations yielding the normalizations by their modular variants as introduced 
by [Arnold 2003] , [Idrees et al. 2011] , and [Bohm et al. 2012a] . These variants 


are either probabilistic or require expensive tests to verify the results at the end. 
In order to reduce the number and complexity of the verification tests, we provide 
a direct modularization for the adjoint ideal algorithm. The approach we propose 
requires only the verification of the final result: We give efficient conditions for 
checking whether the result obtained is indeed the Gorenstein adjoint ideal. 

Our approach relies on the general scheme for modular computations presented 
in | Bohm et al. 2012d|. This scheme is based on error tolerant rational reconstruction 


8.6 


below) and can handle bad prime^ of various types, provided there 


(see Remark 

are only finitely many such primes. Referring to [Bohm et al. 2012d for details, we 


will now outline the main ideas behind the scheme. 

Fix a global monomial ordering > on the monoid of monomials in the variables 
X = {Xo,..., X r }. Consider the polynomial rings R = Q[X] and, given an integer 
N > 2, Rn = (Z/NZ)[X\. If H C R or H C Rn is a Grobner basis, then denote by 
LM (H) := (LM(/) | / € H} its set of leading monomials. 


Iff 


E 


with gcd(a, b) = 1 and gcd(6, N ) = 1, set (|) N := (a + NZ)(b + NZ) 


-1 


Z/IVZ. If / G R is a polynomial such that N is coprime to any denominator of a 
coefficient of /, then its reduction modulo N is the polynomial /jv E Rn obtained by 
mapping each coefficient x of / to xn • If H = {h \,..., ht} C R is a Grobner basis 
such that N is coprime to any denominator in any hi, set Hn = {(/ii)jv, • • •, (^t)iv}- 
If J C R is an ideal, we write 

Jo = J fl Z[X] and Jn = (In \ f E Jo) C Rn, 

and call Jn the reduction of J modulo N. We also write ( R/J)n = Rn/Jn- 
Based on this notation, we fix the following setup for the rest of this section: 

Notation 8.1. Let T C Pq be a curve of degree n. As before, suppose that T is 
integral and non-degenerate. Denote by I(r) the ideal of T in R, and by (7(0) C R 
the reduced Grobner basis of 0(T). If p is a prime such that J(r) p is radical and 


^In our context, a prime p is bad if Algorithm [3] applied to the modulo p values of the input over 
the rationals does not return the reduction of the characteristic zero result. 
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defines an integral, non-degenerate curve in PJ. , then write T p for this curve and 
G(p) C R p for the reduced Grobner basis of 0(r p ). 


Remark 8.2. Given p, the ideal I{ T) p can be found using Grobner bases over Z 
(see |Adams and Loustaunau 1994, Cor. 4.4.5] and [Arnold 2003, Lem. 6.1]). We 


will make use of this in the final verification test. With regard to the other steps of 
our algorithm (in particular, in a randomized version of the algorithm obtained by 
omitting the verification test), we can proceed in the following, more efficient way: 
Let {/i,..., f r } be the reduced Grobner basis of I(T). Reject p if one of the (fi) p is 
not defined (there are only finitely many such primes p). Otherwise, realize I(T) P 
via the equality 


I(T) p = ((h) p ,...,(fr) P )cR p , (8.1) 

which holds true for all but finitely many primes p. These finitely many bad primes 
will not influence the lift if we apply error tolerant rational reconstruction as de¬ 
scribed in Remark 18.61 below. 


Remark 8.3. There are only finitely many primes p for which the desired conditions 
on I(T) P in Notation 8.1 are not satisfied. Since these conditions can be checked 
using polynomial factorization and Grobner bases, we may simply reject such a bad 
prime if we encounter it in our modular algorithm. Hence, we will ignore these bad 
primes in the following discussion. In particular, we will assume that the Grobner 
bases G(p) exists for all primes p. 


The basic idea of the modular adjoint ideal algorithm can then be described 
as follows: First, choose a set of primes V and compute G(p) for each p E V. 
Second, lift the G(p) coefficientwise to a set of polynomials G C R. Provided that 
0(T) P = <25 (r p ) for each p G V, we then expect that G is a Grobner basis which 
coincides with our target Grobner basis G(0). 

The lifting process consists of two steps. First, use Chinese remaindering to lift 
the G(p) C R p to a set of polynomials G(N ) C Rn, with N := Y\ p& -pP • Second, 
compute a set of polynomials G C R by lifting the coefficients occurring in G(N ) 
to rational coefficients. Here, to identify Grobner basis elements corresponding to 
each other, we require that LM(G(p)) = LM(G(g)) for all p,q E V. This leads to 
condition (L2) in the definition below: 

Definition 8.4. With notation as above, a prime p is called lucky if: 

(LI) 0(r) p = 0(r p ) and 
(L2) LM(G(0)) = LM(G(p)). 

Otherwise p is called unlucky. 


Lemma 8.5. All but finitely many primes are lucky. 

Proof. As is clear from the proof of |Bohm et al. 2012d, Lemma 5.5], it is enough to 
show that condition (LI) is true for all but finitely many primes. For this, we may 
assume that both T and do not have any singularities at Xq = 0. Let C be the 
affine part of T. Write A = Q[Xi, ..., X r \/I(C). As shown in [Bohm et al. 2012a], 
(A)p = A p for all but finitely many primes p. So if we write A = fU, with an ideal 
U C A and an element 0 d £ A, and A p = U(p), with U(p) C A p and d p € A p , 
then 

(d p : Up) = (d(p) : U(p)) 

for all but finitely many primes p. Computing an ideal quotient amounts to a 
Grobner basis computation. Hence, as pointed out in Bohm et al. 2012d, Remark 
5.3], 

(d : U) p = (dp : Up) 


for all but finitely many primes p. The result follows, thus, from Propositions |3.9| 
and 14.81 □ 
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When performing our modular algorithm, condition (LI) can only be checked 
a posteriori: We compute G(p) and, thus, ®(r p ) on our way, but ®(r) p is only 
known to us after (7(0) and, thus, ©(T) has been computed. This is not a problem, 
however, since the finitely many primes where ®(T) p ^ ®(r p ) will not influence the 
final result if we apply error tolerant rational reconstruction and the set V is large 
enough: 


Remark 8.6. Let N' and M be integers with gcd(iV / , M ) = 1, let N = N' ■ M, and 
let § G Q with gcd(a, b) = gcd(iV / , b) = 1. Set r\ := (| ) N , G Z/IV'Z, let G Z/MZ 
be arbitrary, and denote by r the image of (r±,r 2 ) under the isomorphism 

Z/7V'Z x Z/MZ -> Z/7VZ. 

Lifting r to a rational number by Gaussian reduction, starting from (ao,£>o) = 
and (ai,&i) = (r, 1), we create the sequence (a,i,bi) obtained by 

(®i +2 5^1+2) — ifliibi) Qi (n *+1 , &i+l), 


with 


Qi = 


Computing this sequence until 


||(aj + i, 6 i + i)|| L > N, and otherwise. By 


((cij, 6j), (dj-j-i, 6j + i)) 

\\{a i+1 ,b i+1 )\\ 2 
(a i+ 2 ,b i+ 2 )\\ > ||(at-n , &i+i)| |, we return false if 


Bohm et al. 2012d 


Lemma 4.3], this 


algorithm will return provided that N is large enough and M <C N'. More 

precisely, we ask that N' > (a 2 + b 2 ) ■ M. 


Definition 8.7. If V is a finite set of primes, set 

N 1 = p and M = p. 

pEV lucky pEV unlucky 

Then V is called sufficiently large if 

N' > ( a 2 + b 2 ) • M 

for all coefficients | of polynomials in G( 0 ) (assume gcd{a,b) = 1 ). 


Lemma 8.8. If V is a sufficiently large set of primes satisfying condition (L2), then 
the reduced Grobner bases G(p), p G V, lift to the reduced Grobner basis (7(0). 


Proof. See [Bohm et al. 2012d, Lemma 5.6]. □ 

From a theoretical point of view, Lemma |8.5| guarantees that a sufficiently large 
set V of primes satisfying condition (L2) exists. From a practical point of view, 
however, (L2) can only be checked a posteriori. Nevertheless, in order to be able to 
identify Grobner basis elements in the lifting process, we have to restrict to a set 
of primes p which all have the same associated set of lead monomials LM(G(p)). 
Hence, taking Lemma 8.5 into account, we proceed along the following lines: First, 
fix an integer t > 1 and choose a set of t primes V at random. Second, compute 
GV = {G(p) | p G V} and use a majority vote with respect to (L2): 


deleteByMajorityVote: Define an equivalence relation on V by setting p ~ 
q :•<=>■ LM(G(p)) = LM(G(g)). Then replace V by the equivalence class of largest 
cardinality jj and change GV accordingly. 

Now, all G(p), p G V, have the same set of leading monomials. Hence, we can apply 
the rational reconstruction algorithm to the coefficients of the Grobner bases in GV. 
If this algorithm returns false at some point, we enlarge the set V by t primes 
not used so far, and repeat the whole process. Otherwise, the lifting yields a set of 

®We have to use a weighted cardinality count: when enlarging V, the total weight of the elements 
already present must be strictly smaller than the total weight of the new elements. Otherwise, 
though highly unlikely in practical terms, it may happen that only unlucky primes are accumulated. 
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polynomials G C R. Furthermore, if V is sufficiently large, all primes in V satisfy 
condition (L2). Since we cannot check, however, whether V is sufficiently large, a 
final verification step is needed. Since this may be expensive, especially if G G(0), 
we first perform a test in positive characteristic: 

pTest: Randomly choose a prime p (/ V which does not divide the numerator 
or denominator of any coefficient occurring in a polynomial in G. Return true if 
G p = G(p), and false otherwise. 

If pTest returns false, then V is not sufficiently large (or the extra prime chosen 
in pTest is bad). In this case, we enlarge V as above and repeat the process. If 
pTest returns true, however, then most likely G = G(0). In this case, we verify 
the result over the rationals as described below. If the verification fails, we again 
enlarge V and repeat the process. 

We now discuss the verification. We write I = {G)r for the lifted modular result 
and © = ®(r) C R for the correct result. After checking that G is indeed a Grobner 
basis and / is saturated (henceforth, this will be assumed), we apply the following 
results. 


Lemma 8.9. With notation as above, the ideal I is equal to the Gorenstein adjoint 
ideal & of T iff 


(1) m g 1, 

(2) deg A(/) = deg/+ d(r), and 

(3) deg I = deg <5. 


Proof. If 7 = 0, then I satisfies (1), (2), and (3). Conversely, by Lemma 3.13 


conditions (1) and (2) imply that I is an adjoint ideal of T. In this case, since © is 
the largest such ideal, we have / C ©. But then I = & by (3). □ 


It is clear how to check condition (1). In what follows, we describe a method 
for checking (2) which, in particular, provides a way of finding deg A(0). This will 
allow us to check (3) via the formula deg© = degA(©) — <5(r). 

If k is any field, and A is any reduced Noetherian /-algebra, the delta invariant 
of A is defined to be 

5 k (A) = dim*. A/A. 


Proposition 8.10. Let B be a ring, and let A be a B-algebra with the following 
properties: 

(1) (B, m) is a normal local ring with perfect residue class field k. 

(2) B -A B is fiat, and for all p £ Spec(/>) such that pB B, the ring B ® 
Bp/pB v is geometrically normal. 

(3) A is a formally equidimensional Nagata ring. 

(4) A is a flat B-algebra, md is contained in every maximal ideal of A, A/mA 
is reduced, and 8 k {A/mA) < 00 . 

(5) A/A is a finite B-module. 

(6) The unique map A/mA —> A/mA factorizing the normalization map A/mA — > 
A/mA as 

A/mA -a A/xnA -a A/mA 

is injective. 


Then 


Sq (b){A ®b Q (B)) < S k (A/mA). 


Proof. See |Lipman 2006, Prop. 2.2.1(i)] for the factorization in (6) and Lipman 


2006, Prop. 3.3] for the proof of the proposition. 


□ 


Corollary 8.11. In the setting of N otation\8 . 1\ given a prime p, we have 


m < s(t p ). 
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Proof. Let X' = {Xi,...,X r }. We may assume that T has no singularities at 
Xq = 0. As before, let C be the affine part of L. Then J := 1(C )o C Z[X'\ is a 
prime ideal of height to —1, (p, J) is a prime ideal, and Jn Z = (0). The claim follows 
by applying Proposition |8.10 to (B, m) = (Z^,(p)) and A = 


since, then, A® B Q(33) = Q[X']/3(C) and A/mA = F p [A']//(C') p , and conditions 
(1) through ( 6 ) of the proposition are satisfied. Indeed, this is clear for (1), and (2) 
holds since B is excellent. Moreover, we have (3) since A is of finite type over B and 
J7j^[X'] is a prime ideal. Condition (4) follows since A is a torsion free 13-module, 
(p, J) is a prime ideal, and Spec(A/mA) is a curve. We obtain (5) since A/Ca is 
a finite 33-module and A/Ca is a finite A/C^-module. Condition ( 6 ) follows from 
Lemma 8.13 below which gives us a canonical map 

a mod (p, J) 


A —> A/mA, 


a = = >-a 
b 


b mod (p, J) ’ 


where a, b are the images of a, b E [. X'] in A, and b (p, J). Since a = | is in 

the kernel of this map iff a E {p, J), we get an injective map A/mA —> A/mA which 
factors the normalization map as desired. □ 


Before deriving Lemma 8.13, we illustrate condition ( 6 ) by an example. 
Example 8.12. Let (33,m) = (Z ( 3 ) ,<3» and A = Z ( 3 ) [X, Y\/ (X 3 + Y 3 + T 5 ). 
Then A/mA = ^1, |, ^ and A = ^1, |, ^ We compute <Sq(A<8 > b Q) = 

3 and Sp (A/mA) = 4, and find that 


A/mA = = A/mA. 

y y / A/mA \ 2/ r / A/mA 


Lemma 8.13. Wif/i t/ie notation of the proof of Corollary 8.11. for any a E A f/iere 
exist a, 6 E Z[W] with b ^ (p, J) such that a = |. 


Proof. For a e A C Q(A) = Q(Z[X']/J), there are a,b £ r L\X'\ with b ^ J and 
a = j , and there are ao,..., a m _i E Z[X 7 ] and d E Z with p f d such that 


a m + am-1 m °d d ^ m _i 


an mod J 

... + JL -j— =0, 


that is, d ■ a m + a m _i • ba m ~ l + ... + ao ■ b m E J. 

If b E (p, J), then d • a m E (p, J), hence, since J is radical, a E (p, J). Then 
a = pai + ci and b = pb\ + di with ai, 6 i E Z[X 7 ] and ci, di E J. If 6 i E (p, J), we 
can iterate the process. Inductively, we obtain a s ,b s E Z[A' 7 ] and c s ,d s E J with 
a = p s a s + c s and b = p s b s + d s . If b s E (p, J) for all s, then 6 E (p s , J) = J, a 
contradiction. Otherwise there is an s with b s (p, J). Then 


amodJ p s a s mod J a s modJ 
b mod J p s b s mod J b s mod J 


□ 


In the following, we write again it : T — > T for the normalization map, and denote 
by M the vanishing ideal of Sing(T) in R. Consider a homogeneous polynomial 
g E 3 = {G) r not contained in 3(T), and let m be its degree. Let div(p) be the 
divisor cut out by 7 r *g on T, let D(g) = div(p) — A(3) be the corresponding divisor 
in |m33 — A(3)|, and let d(g) = deg D(g). Furthermore, write d(g) for the degree of 
the part of D{g) away from Sing(r). Then d{g) < d(g), and d(g) can be computed 
as 

d(g) = deg ((3(T) + (g)) : M°°) 
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provided that / : M°° = (1), what we will henceforth assume (in Algorithm [ 6 ] below, 
if this condition is not fulfilled, we enlarge our set of primes). 

Theorem 8.14. Let I = {G)r be as above, and let p be a prime number. Suppose: 

( 1 ) LM(/(r p )) = LM(/(r)), 

(2) G(p) is a Grobner basis of an adjoint ideal of T p , 

(3) G p = G(p), 

(4) d(g p ) = deg(r) • m - deg((G(p)) flp ) - 5(T), and 

(5) m is large enough to ensure that \mH — A(/)| is nonspecial. 

Then 

deg A (!) = deg(r) • m - d(g p ). 

Furthermore, degA(I) = degA(/ p ), and I is an adjoint ideal ofT. 

Remark 8.15. To apply the theorem in the setup above, note: Condition (1) can 
easily be tested. Furthermore, ( 2 ) and (3) are satisfied by the construction of G. 
Since we know how to compute <5(r), condition (4) can be tested. With respect to 
(5), we will comment on how to choose m in Lemma 8.17 below. 

Proof of the theorem. By (1), deg(r p ) = deg(T) and p a (T P ) = p a (r). First note, 
that by (3) 

I P = ( G)r p = (G(p)) Rp 

and, as G is assumed to be a Grobner basis, 

deg«G) fl ) = deg((G(p)) fl ). (8.2) 


By Corollary 8.11 we have 5(F) < <5(r p ). Hence 

d{g P ) < d(g p ) = deg(r p ) • m - deg A(/ p ) 

= deg(r) • m - deg (Ip) - 5(r p ) 

< deg(r) • m - deg(/ p ) - <5(r) 

using that by (2) the ideal I p is an adjoint ideal of T p . By (4) the chain of inequalities 
is an equality, hence 

d{g P ) = d(g p ) = deg(r p ) • m - deg A(/ p ) 

and 

*(F) = S(T p ). 

By (8.2) and Lemma |3.13 this implies that 

deg A {I p ) = deg(Jp) +5(T p ) = deg (!) + <5(r) > degA(I), (8.3) 

or equivalently 

d{g p ) < d(g). 

To prove equality, we consider the closed subscheme 

X = V(/(r)o) cPg A SpecZ 

with projection n and fibers X q = n~ 1 ((q}) = X xg pec z Spec K((q)). So over the 
generic point (0) G SpecZ the fiber is Xq = F and over (p) it is X p = F p . By (1) 
the Hilbert polynomials of T and T p are equal, hence there is a Zariski open subset 
V C SpecZ with (p) G V such that the Hilbert polynomial is constant on V. So 
Try : Xy = 7 r _ 1 (H) —» V is a flat family (see |Hartshorne 1977, Ch. Ill, Thm. 9.9]). 

Since <5(r p ) = 5(T), the 5-constant criterion for simultaneous normalization (see 
|Lipman 2006]) implies that there is a Zariski open subset U C V C SpecZ with 
{p) G U such that nu is equinormalizable. That is, there is a finite map v : Z -A Xjj 
such that 7f := 7 tjj o v is flat with nonempty geometrically normal fibers, and for 
each (q) G U the induced map on the fibers v q : X q = 7 f - 1 ((g)) —> 7 r^ 1 ((g)) = X q is 
a normalization map. 
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Since, by construction, the family of sheaves defined by Iq is flat over U and 


U contains both (0) and (p), the semicontinuity theorem (see, for example, Liu 


2002, Ch. 5, Thm. 3.20]) implies that the dimensions of the linear series induced 
by I on T and by I p on T p satisfy 

h° (f p , 0 Tp {m • H p - A (Ip))) > h° (f, O r (m • H - A (/))) . 

Hence by (5), Riemann-Roch, and <5(T p ) = <5(T) it follows that the degrees of the 
linear series satisfy d(g p ) > d(g), so we obtain the second equality in 

d(g P ) = d{g p ) = d(g) 

(having shown the first already above). The second equality also translates into 
deg A (I p ) = deg A(I) which, by ( |8.3[ ), implies that I is an adjoint ideal. Moreover, 

deg(T) - m- deg A(J) = deg(T p ) • m - deg A (I p ) = d(g p ). 

□ 


Remark 8.16. Suppose now, in addition to the previous assumptions, that I p is the 
Gorenstein adjoint ideal of T p . Since / is an adjoint ideal of T, we have I C © which 
implies deg I > deg0, hence 

degA(C5) = deg(6) + <5(r) < deg(I) + 5(T) = degA(I) = degA(/ p ). (8.4) 
Moreover, by semicontinuity 


dim |m ■ H p — A(I P )| > dim |m ■ H — A(©)| 
for m large enough, so by Riemann-Roch and £(T p ) = <5(T) we have 

deg A(/ p ) < deg A(0). 


Hence (8.4) is an equality and implies 

deg / = deg C5, 


that is, I is the Gorenstein adjoint ideal of T. 


In order to expect condition (4) to be satisfied for randomly chosen g and p, the 
degree m has to be chosen large enough such that d(g) = deg(T) • m — degA(©) 
for a generic g € © m (taking into account that d(g) = d(g p ), and <5(T) = 5(T p ) and 
Ip = hence deg A(/ p ) = degA(6) holds true for all but finitely many primes 

p ). The following lemma specifies an appropriate bound for m, which will also be 
sufficient to obtain (5). 

Lemma 8.17. Consider an integer m such that Pr(m) — 1 > p a (T) and suppose 
that g G © m is generic. Then 

d(g) = deg(T) • m — deg A(0) 

Furthermore, \mH — A(©)| is nonspecial. 


Proof. By assumption and since Pp(m) = (degT) ■ m — p a (T) + 1, we have 

deg(r) • m > 2p a (T). 

By Corollary 3. Til] we obtain degA(C5) < 2d(T). Hence, it follows that 
deg(T) • m — deg A(0) > deg(r) • m — 2<5(r) 

= deg(T) • m - 2p a (T) + 2p(T) > 2p(T). 

This implies that \mH — A(0)| is base-point free (see |Hartshorne 1977, Ch. IV, 
Cor. 3.2]), hence, since g is generic, we have d(g) = d(g). By reason of its degree, 
the linear series is also nonspecial (see Hartshorne 1977, Ch. IV, Ex. 1.3.4]). □ 
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Remark 8.18. For a plane curve T of degree n the condition Pp(m) — 1 > p a { r) is 
equivalent to n ■ m > (n — l)(n — 2), which is satisfied for m > n — 2. 

We summarize our approach in Algorithm [6j 


Algorithm 6 Modular adjoint ideal 


Input: A curve T C P r satisfying the conditions of Notation 8.1 
Output: The Gorenstein adjoint ideal 0(T). 

1: choose an integer t > 1 

V = QV = 0 

loop 

choose a list Q of t random primes not used so far 

for all p G Q do 


9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 


if T p is irreducible, non-degenerate, and LM(/(T)) = LM(/(T P )) then 
compute the reduced Grobner basis G p of ®(T p ) C R p (via Alg. [3]) 

V = V U {p}, QV = QV U {G p } 

(QV,V) = deleteByMajorityVote^R, V) 

lift {QV, V) to a set of polynomials G C R via the Chinese remainder theorem 
and Gaussian reduction 

if the lifting succeeds and pTest(/(T), G, V) then 

if G is a Grobner basis and (G) is saturated and (G) : M°° = (1) then 
choose m such that Pp(m) — 1 > p a (T) 
choose g G (G) m at random 
choose a prime p G V 

compute d{g p ) = deg ((/(r p ) + (g p )) : M£°) 


compute <5(T) by applying Remark 7.3 


if d(g p ) = deg(r) • m — deg ( G(p )} — <5(r) then 
return (G) 


Remark 8.19. In Algorithm [6j the different G{p) can be computed in parallel. The 
individual computations can be parallelized by partitioning the singular loci. 

Remark 8.20. The most expensive step of the verification is the computation of 
<5(r). If we skip the verification, the algorithm will become probabilistic, that is, 
the output is the Gorenstein adjoint ideal only with high probability. This usually 
accelerates the algorithm considerably and gives us, in particular, a fast probabilistic 
way to compute both the geometric genus p(r) and deg A(0) = dim^ 


c] 


9. Timings 


The algorithms for adjoint ideals presented in this paper are implemented in the 
SINGULAR library adjointideal.lib (see Bohrn et al. 2015b|). They make use 
of the normalization algorithm of Section [2] either in its local or local to global 
variant, as appropriate. These variants, in turn, are part of the Singular library 

locnormal.lib (see 


Bohrn et al. 2012b 


In this section, we compare the performance of the different algorithms. Specifi¬ 
cally, we consider 


LA 

IQ 

locIQ 

locIQP2 


modLocIQ 


Mnuk’s global linear algebra approach (Algorithm [lj) , 

the global ideal quotient approach (Algorithmpl), 

the local ideal quotient approach (Algorithm |3]using Algorithm |4j) , 

the local ideal quotient approach for plane curves with the 

improvements of Section Pd concerning ordinary multiple 

points and singularities of type ADE, and 

the modular local ideal quotient strategy (Algorithm [bj) . 














28 


J. BOHM, W. DECKER, S. LAPLAGNE, AND G. PFISTER 


For the modular approach, we do not make use of a local analysis of the singular 
locus except for computing the invariants needed in the verification step. 

To quantify the improvement in computation time obtained by omitting the veri¬ 
fication step in the modular approach, we give timings for the resulting, now proba¬ 
bilistic, version of Algorithm^ (denoted by modLocIQ J in the tables). In all examples 
computed so far, the result of the probabilistic algorithm is indeed correct. 

To quantify the contributions of the different normalization algorithms and to 
provide a lower bound for any adjoint ideal algorithm using them, we also specify 
the following computation times: normalization in Singular via the local to global 
approach outlined in Section [ 2 ] (denoted by locNormal); and finding an integral basis 
in Maple via the algorithm of van Hoeij (denoted by Maple-IB). Once being fully 
implemented in Singular, we expect further improvements of the performance by 
computing the local contribution or just an integral basis of the local ring by the 
algorithm discussed in Bohm et al. 2015a]. Since this algorithm and van Hoeij’s 


algorithm rely on Puiseux series, they work in characteristic zero only. 

All timings are in seconds on an AMD Opteron 6174 machine with 48 cores, 
2.2GHz, and 128GB of RAM running a Linux operating system. A dash indicates 
that the computation did not finish within 10000 seconds. The timings for parallel 
computations are marked by the symbol * and the maximum number of cores used 
in parallel is indicated in brackets. 


Remark 9.1. All examples are defined over the field of rationals. For locIQ*, the 
number of cores used corresponds to the number of components of the decomposition 
of the singular locus over Q. For modLocIQ*, the number of cores used in a given 
iteration of the algorithm is obtained by summing up the number of components 
modulo p over all primes p E Q chosen in Step 4 of Algorithm [6j 

To show the power of the modular algorithm, we give simulated parallel timings 
even if the number of processes exceeds the number of cores available on our machine 
(which is a valid approach since the algorithm has basically zero communication 
overhead). For the single-core timings of modLocIQ, we indicate in square brackets 
the number of primes used by the algorithm. 

Now we turn to explicit examples. First we consider rational plane curves de¬ 
fined by a random parametrization of degree n. These curves have ("g 1 ) ordi¬ 
nary double points. Their defining equations fi, n were generated by the function 
randomRatCurve from the SINGULAR library paraplanecurves. lib (see [Bohm et 
al. 20 12c| ), using the random seed 1 and a random parametrization with coefficients 
of bitlength 15. 



/l,5 

/l,6 

fl,7 

deg 

5 

6 

7 

locNormal 

2.1 

56 

— 

Maple-IB 

5.1 

47 

318 

LA 

98 

4400 

— 

IQ 

2.1 

56 

— 

locIQ 

1.3 

54 

3800 

locIQ* 

1.3 (1) 

54 (1) 

3800 (1) 

locIQP2 

.18 

1.2 

49 

locIQP2* 

.18 (1) 

1.2 (1) 

49 (1) 

modLocIQ 

6.4 [33] 

19 [53] 

150 [75] 

modLocIQ ’ 

6.2 [33] 

18 [53] 

104 [75] 

modLocIQ* 

.36 (74) 

1.6 (153) 

51 (230) 

modLocIQ 1 * 

.21 (74) 

0.48 (153) 

5.2 (230) 


We observe that the detection of special types of singularities is fast and yields the 
best performance among the non-probabilistic algorithms. 
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To compare the algorithms at a single singularity, we consider plane curves with 
exactly one A n respectively D n singularity at the origin of the affine chart {Z / 0} 
(ignoring singularities at infinity). For the modular approach, we omit verification 
since this step relies on global properties of the curve. 

The curves with affine equation f'z.n.d = Y 2 + X n+1 + Y d , n > 1, d > 3, have 
precisely one singularity of type A n at the origin: 



f2,5,10 

/ 2 , 5,100 

f2,5,500 

/ 2 , 50,100 

/ 2 , 50,500 

/ 2 , 400,500 

deg 

10 

100 

500 

100 

500 

500 

locNormal 

.12 

.12 

.12 

.51 

.51 

3.6 

Maple-IB 

.08 

1.5 

96 

4.7 

150 

630 

LA 

.18 

140 

- 

150 

- 

- 

IQ 

.12 

.12 

.12 

.51 

.51 

3.6 

modLocIQ’ 

.20 [2] 

.22 [2] 

.96 [2] 

1.1 [2] 

2.0 [2] 

11 [2] 

modLocIQ’* 

.10 (2) 

.13 (2) 

.48 (2) 

.54 (2) 

1.2 (2) 

5.8 (2) 


The curves with affine equation f = X(X n 1 + Y 2 ) + Y d , n > 2, d > 3, have 
exactly one singularity of type D n at the origin: 



1 * 3 , 5,10 

1 * 3 , 5,100 

1 * 3 , 5,500 

/ 3 , 50,100 

/ 3 , 50,500 

fa, 400,500 

deg 

10 

100 

500 

100 

50 

500 

locNormal 

.15 

.15 

.15 

.67 

.67 

4.9 

Maple-IB 

.05 

1.7 

100 

34 

1830 

- 

LA 

.20 

140 

- 

140 

- 

- 

IQ 

.15 

.15 

.15 

.67 

.67 

5.0 

modLocIQ ’ 

.22 [2] 

.23 [2] 

.23 [2] 

1.5 [2] 

1.5 [2] 

24 [2] 

modLocIQ’* 

.09 (2) 

.10 (2) 

.10 (2) 

.74 (2) 

.77 (2) 

17 (2) 


In both examples, the best strategy is IQ since we consider only one singularity 
and since no coefficients of large bitlength occur. 

The plane curves with defining equations 


f A n = ( X n+1 + Y n+1 + Z n+1 ) 2 - 4 ( x n+1 Y n+1 + y n+1 Z n+1 + z n+1 X n+1 ) 

were given in |Hirano 1992) and have 3 (n + 1) singularities of type A n if n is even. 
To ensure that all singularities of the curves are in the affine chart {Z ^ 0}, we 
substitute Z = 2X — 3Y + 1. 



/4,4 

.A,6 

OO 

< 

deg 

10 

14 

18 

locNormal 

1.6 

- 

- 

Maple-IB 

2.2 

14 

70 

LA 

89 

- 

- 

IQ 

2.5 

- 

- 

locIQ 

.96 

- 

- 

locIQ* 

.36 (6) 

- 

- 

iociQP2 

1.0 

- 

- 

iociqp2* 

.38 (6) 

— 

— 

modLocIQ 

3.7 [3] 

23 [4] 

190 [4] 

modLocIQ’ 

3.3 [3] 

20 [4] 

170 [4] 

modLocIQ* 

.63 (27) 

4.4 (48) 

50 (48) 

modLocIQ’* 

.38 (27) 

2.2 (48) 

30 (48) 


To conclude this section, we present examples of curves in higher-dimensional 
projective space. As above, we first consider curves with only one singularity in a 
given affine chart: let L n be the ideal of the image of 
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Second, denote by I n the ideal of the image in P 5 under the degree-2 Veronese 
embedding of the curve {/ 4 , n = 0}. The resulting timings are: 



T 25 

£50 

h 

h 

deg 

25 

50 

20 

28 

locNormal 

3.9 

84 

21 

- 

IQ 

3.9 

84 

30 

- 

locIQ 

3.9 

84 

18 

— 

locIQ* 

3.9 (1) 

84 (1) 

7.5 (6) 

- 

modLocIQ’ 

6.5 [2] 

220 [2] 

74 [5] 

2600 [5] 

modLocIQ’* 

3.3 (2) 

140 (2) 

4.0 (45) 

59 (69) 


To summarize, we observe that the ideal quotient approach is faster than the 
linear algebra one. To some extent, this is due to the lack of efficiency of the 
rational function arithmetic in Singular. The local strategy is faster than the 
global one if there is more than one component in the decomposition of the singular 
locus over Q. In addition, the local algorithm can be run in parallel and is, then, 
even faster. In most examples, especially when the coefficients have large bitlength, 
the fastest approach is the modular local strategy, which parallelizes in a two-fold 
way, by localization and modularization. In contrast to other modular algorithms 
(such as modular normalization), the verification step is usually very fast. 

Acknowledgements. We would like to thank Christoph Lossen, Thomas Markwig, 
Mathias Schulze, and Frank Seelisch for helpful discussions. 
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